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AN APPROXIMATE SOLUTION OF ADDITIVE-DRAG COEFFICIENT 
AND MASS- FLOW RATIO FOR INLETS UTILIZING RIGHT 
CIRCULAR CONES AT ZERO ANGLE OF ATTACK 

By Vincent R. Mascitti 
Langley Research Center 

SUMMARY 

This paper presents an approximate solution of additive-drag coefficient and mass- 
flow ratio for inlets utilizing right circular cones at zero angle of attack. The results 
are in good agreement with the exact solution, especially at high Mach numbers and large 
cone half-angles. The present method uses the incompressible approximation of the 
conical flow field. This conical-flow-field approximation facilitates analytic integration 
of the pressures acting on the entering streamline. In addition, the flow-field properties 
obtained by the approximation are in good agreement with the results of the exact solu- 
tion and other approximate solutions. The largest deviation between the exact solution 
and the present solution of additive-drag coefficient was AC-p a = 0.039. 

INTRODUCTION 

The calculation of additive-drag coefficient and mass-flow ratio for axisymmetric 
inlets is complex and time-consuming. The exact solution may be obtained by numerical 
integration of the pressures acting along the entering streamline. This procedure was 
used to obtain the results presented in reference 1. Using graphical or tabulated data 
for automated computation of installed performance of air-breathing engines requires an 
elaborate system of fitted curves over a wide range of conditions. This approach may 
limit computer storage, generality, and, most important, accuracy. 

This paper presents an approximate method by which additive -drag coefficient and 
mass-flow ratio can be calculated directly. Direct calculation has been accomplished by 
using an approximate solution to represent the conical flow field. Use of this approxima- 
tion permits the force on the curved streamline to be determined by analytic integration. 
The approximate expressions for additive-drag coefficient and mass-flow ratio for inlets 
utilizing right circular cones at zero angle of attack have been compared with the exact 
solution of reference 1. 



SYMBOLS 


a 

nondimensional speed of sound, a/Vj 

a 

speed of sound 

A 

cross-sectional area 

c D,a 

D a 

additive -drag coefficient, — ~ 

qAL 

c p 

pressure coefficient 

D a 

additive drag 

M 

Mach number, v/a 

P 

pressure 

q 

dynamic pressure, | p^M^ 


u nondimensional velocity along conical ray line in spherical coordinate 

system, u/vjj 

v nondimensional velocity normal to conical ray line in spherical coordinate 

system, vjvi 

u,v velocity components, dimensional 

V nondimensional resultant velocity at any ray line, v/vi 

V resultant velocity, dimensional 

Vi limiting velocity due to adiabatic expansion into a vacuum 

w mass flow 

w/Woq mass -flow ratio 

y ratio of specific heat at constant pressure to specific heat at constant volume 
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e 


conical ray angle measured from cone axis 


Subscripts: 

s cone surface 

L location of inlet lip 

t stagnation conditions 

w conditions back of shock wave 

°° free-stream condition 

METHOD OF SOLUTION 
Conical- Flow Approximation 

The purpose of this section is to derive a conical-flow-field approximation which 
permits the analytic integration of the pressures acting on the entering streamline. Also, 
the conical-flow-field approximation should provide good agreement with the exact solu- 
tion over a wide range of conditions. 

Two conical-flow approximations which satisfy the latter condition are presented 
in references 2 and 3. The method of reference 2 offers little advantage over the exact 
solution because it requires an iterative or semigraphical solution. In contrast, the 
method of reference 3 is a closed-form solution. Unfortunately, attempts to integrate 
the pressures acting on the entering streamline failed. 

A conical -flow approximation has been found which satisfies the two requirements 
stated above. The differential equation which corresponds to this solution was first 
recognized in reference 4 as the incompressible solution to the supersonic conical flow 
field. This solution is more appropriately referred to as a constant-density solution as 
pointed out in reference 5. In reference 6, the solution to this equation has been applied 
to the hypersonic speed regime where the constant -density approximation is valid. How- 
ever, it is not widely recognized that this solution gives good results even at the lower 
Mach numbers. The "Results and Discussion" section indicates that the results of the 
constant-density solution are comparable to those obtained by the approximation of 
reference 3. 

The derivations of the basic equations for the conical -flow problem are presented 
in reference 4 and, therefore, will not be repeated herein. The second-order, nonlinear 
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differential equation representing the conical flow field in the spherical coordinate sys- 
tem (fig. 1) is 

d 2 u + u _ a 2 (u + v cot g) /-j 

d0 2 v 2 - a2 


where 

and 


In these equations, the velocities had been nondimensionalized by dividing them by the 
limiting velocity attainable by adiabatic expansion into a vacuum. Rearranging equa- 
tion (1) gives 

O 

d^u u + v cot 6 
d0 2 v 2 

a 2 

The equation is nonlinear because of the quantity v 2 /a 2 . The assumption of this paper 
v 2 

is that — « 1 (constant-density solution). 
a 2 

The differential equation reduces to 

2 

iJi + cot 6 ^ + 2u = 0 (2) 

d$2 d6 


v = 


du 

d0 


a 2 = 21_I_i(i - u 2 - v 2 ) 
2 


Equation (2) is linear and has as one solution u = cos 6. With this solution, the second 
solution can be found by quadratures, as follows: 


where and C 2 
conditions. 


u = cos 9<C 1 


sec 6 + 


l°g e (tai|) 


+ c. 


are constants of integration to be evaluated by the boundary 


(3) 


One boundary condition is at the conical surface when 6 = 6 S and is as follows: 


v = 


du 

d6 


0 


This condition results in the relation 


C 2 = Cj 


cos 8c 


sin 2 0 s 


ioge^tan x 


(4) 
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Substituting equation (4) into equation (3) gives 


/ o\ cos 8 S ( 6, 

u = Ci cos 8 sec 6 + log tan - + — log tan 

e \ w sin 2 0s e V 2 


The second boundary condition is behind the shock wave when 6 = 0 W and is as 
follows: 

u w = cos 0 W 

From reference 7, 



Therefore, 



Substituting equation (6) into equation (5) gives 



Substituting equation (7) into equation (5) gives 



where 


( 7 ) 


( 8 ) 


K = 



5 



In addition, 


v = = K Vcot 6 - sin 6\ 


d6 




cos 0 S 

sin^^g 



and 




(9) 


( 10 ) 


Equations (8), (9), and (10) are the expressions for the flow velocities in terms of the cone 
half-angle, shock wave angle, free-stream Mach number, and ray angle. 

It should be pointed out that the Rankine-Hugoniot conditions at the shock have not 
been satisfied. The second boundary condition merely specifies that the tangential veloc- 
ity component is constant across the shock. The reason for this approach is to provide 
two methods for determining the shock wave angle: 

(1) To impose the Rankine-Hugoniot condition and develop an expression for shock 

wave angle as a function of free-stream Mach number and cone half-angle 

(2) To use the exact shock wave angle in tabular form to provide more precise 

results 

With equations (8), (9), and (10) in their present form, either method may be used. 

The Rankine-Hugoniot equations can be expressed as: 


tan 6 


w 


- T ~ 1 u w ~ 1 
y + 1 u w v w 


(ID 


Equations (8) and (9) evaluated at the shock wave become 


u 


w 


1 + 


cos Oyj 


y - 1 


M 0 


( 12 ) 


6 



r 


v w = 


cot 6 W - sin 6 


w 


1 + 


Xll 

y " 1 M J 


cos 0 S 


sin^ 6 S 


+ log ( 


/tan — \ 
2 




W/ 


cos 0 S ' tan 

sec + — — + log e . 


y w> 


sin^0 s 




‘an 


(13) 


Substituting equations (12) and (13) into equation (11) and reducing gives 


1 


y + 1 

f 

1 

''N 

2 

< 

1 

/, 0 W \ 
COS 6L / tan 9 \ 



1 + cos 0 W 

— T^ + logJ f- 

sin 2 0s Itan ^ / 



V. 

\ 2 / J 

J 


(14) 


Unfortunately, the shock angle cannot be calculated directly when cone half-angle 
and free-stream Mach number are specified. Since an iteration scheme is required to 
solve equation (12) when 6 S and M m are specified, using the solution of the approxi- 
mate shock angle offers little advantage over using a tabulation of the exact shock angle. 


Additive-Drag Coefficient and Mass-Flow Ratio 

By using the expressions for the velocity components in the conical flow field, the 
equations for additive-drag coefficient and mass -flow ratio can be derived. The detailed 
steps of the derivation are presented in the appendix. The resulting equation for the 
additive -drag coefficient is 


where 


w r 

4K 

P’ w 

(r - 

Wj 

c - 

2 

?Lsl(i 

Op - 

yMoo 2 

Poo ' 


4K Pt.w v(l - V 2 ) 7 " 1 / tan 2 


log 


e 0. 

\tan 


w 


y 


(15) 
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and, from reference 7, 


v 

y + 1 

1 

Xh 

Vw 

[ (r + l)M oo 2 sin 2 0 w [(y - ljM^ 2 + 2j|‘ 

Poo 

_ 2 rM 0O 2 sin 2 e w - (y - 

' !)J 

1 

1 i 

2 (y - l)M oo 2 sin 2 0 w + 2 J 


The resulting equation for mass -flow ratio is 


1 

_w_ _ _v_ sin e w/l - V 2 ^ 7 " 1 
Woo v w sin e U _ v w 2y 


(16) 


From equations (15) and (16), the two obvious end results are satisfied: 


When 6 = 0 W 


When 0 - 6 S 


w 

Woo 


= 1 


C D,a = 0 


v = 0 


-2L=0 

Woo 


C D,a “ C P,s 


RESULTS AND DISCUSSION 
Conical- Flow Approximation 

The purpose of this section is to provide some confidence in the use of the constant- 
density conical-flow solution over the Mach number range. Evidenced by reference 6, 
the constant-density solution has been successfully applied at high Mach numbers. How- 
ever, no comparison of the results of the approximate methods with those of the exact 
solution was made at lower Mach numbers. Figure 2 presents the variation of the pres- 
sure coefficient with conical ray angle for four cones at zero angle of attack and at a free- 
stream Mach number of 2. The solid line corresponds to the exact solution of refer- 
ence 8. The short dashed line corresponds to the approximate solution using equation (10) 
and the exact value of shock angle. The long dashed line corresponds to the approximate 
solution of reference 3. The results of the present solution agree reasonably well with 
the results of the exact solution and those obtained by the approximation of reference 3. 

Additive-Drag Coefficient and Mass- Flow Ratio 

Additive-drag coefficient and mass-flow ratio, equations (15) and (16), have been 
computed for an inlet cowl placed in the conical flow field at zero angle of attack. The 
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results are presented for cone half-angles of 5°, 10°, 20°, and 30° in figures 3, 4, 5, 
and 6, respectively. Additive-drag coefficient and mass-flow ratio are plotted against 
ray angle for three Mach numbers: sonic flow on the cone surface, Moo = 2, and 
Moo = 6. The solid lines correspond to the exact solution obtained by numerical integra- 
tion in reference 1. The dashed lines were computed from the approximate expres- 
sions (15) and (16). Additive-drag coefficient is in good agreement with the exact solu- 
tion, especially at the higher Mach numbers and cone half-angles. The fact that the 
agreement in mass-flow ratio is very good over the full range of conditions indicates 
that the conical-flow approximation provides a good approximation of the entering 
streamline. 

Figures 3 to 6 show that the maximum deviation in additive-drag coefficient occurs 
when the cowl lip is placed on the cone surface. Figure 7 shows the maximum deviation 
in additive-drag coefficient between the exact and approximate solutions. The difference 
between the exact and approximate additive-drag coefficient is plotted against Mach num- 
ber for four cone half-angles. The largest deviation shown is ACj) a = 0.039 at a Mach 
number of 1.12 for 0 S = 10°. The importance of this error in terms of inlet and aircraft 
performance will vary with aircraft configuration and is beyond the scope of this paper. 

CONCLUSION 

An approximate method for calculating additive-drag coefficient and mass-flow 
ratio of inlets utilizing right circular cones at zero angle of attack has been presented. 
The solution uses the constant -density conical-flow-field equations. These equations 
are shown to provide good agreement with the exact cone solution over a wide range of 
variables. Results of the approximate expressions for additive -drag coefficient and 
mass-flow ratio are in good agreement with the results of the exact solution and other 
approximate solutions. The largest deviation between the exact solution and the present 
solution of additive-drag coefficient was ACj^ a = 0.039. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., September 2, 1969. 
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APPENDIX 


DERIVATION OF EQUATIONS FOR ADDITIVE -DRAG COEFFICIENT 


AND MASS- FLOW RATIO 


Additive drag is defined as 


c e 

D a = (p - p^dA 

A 


The additive -drag coefficient is given by 


Substituting equivalent expressions for the parameters into equation (A2) and rearranging 
yields 


’ yM^AL 


£(i 


Because 


p _ p F t,w 

Poo Pt.w Poo 


-£_= (1 - v 2 ) 7-1 


'D,a = “A— ^ J. ^ (l - V 2 ) r_1 dA - (A L - A w ) 

yMoo^AL 


The mass-flow ratio captured by a cowl lip at position L is defined as follows: 


(l - v 2 ) 7 ' 1 


r °° C sin 6 
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APPENDIX — Continued 


where 


1 - V, 


2V-1. 


sin 6x 


From reference 1, 


Substituting the equivalent expression for w/w^ into equation (A5) yields 


t rC sin 


(1 - V 2 ) 


1 

2 V '- 1 


Then 


dA = cos 6 d6 - sin 0 dv + 2 sin 0V dV 


(l - V 2 ) 


2 ) y ‘ 1 v2 


(r - i)(i - v 2 ) r " v 


Applying the expression for v from equation (9) gives 

v cos 6 d6 - sin 6 dv = d0 

sin 6 

Applying the expressions for v and V from equations (9) and (10) gives 

2 sin 6V dV = _ 2K d& 
v sin 0 

Substituting the equivalent expressions (A6), (A7), and (A8) into equation (A4) gives 


'D,a ~ 


YMj Ar p 


2K(l - V 2 ) 2K 

v 2 sin 0 (7 ~ l)sin 


d0 - (A L - A w ) 


The procedure for solving the first integral is as follows: 


:(i - v 2 ) 


v^sin 


1 - K 2 csc 2 ( 


2 cos 


Jcos 0„ / tan |\ 

-K 2 — ^ + log J -2- 

sin 2 0 s e s I 

s \tan 

= — 4— 2K csc 2 0 d0 

,.-."12 


cos 0g / 

^ + 1 °M 
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APPENDIX - Continued 


R = K 


cos 0 
sin 2 0 


cos 0 S i ( LcU1 2 

S‘" 2 *s + 0 S L^ 


dR = -2K csc 3 0 d0 


2K csc 6 e d6 


2 ) cos 9 


cos e s / tan 2 
sin2e s ge itan ^ 


^ dR = 1 0 

*0 W R 2 R 0 w 


cos 0 
sin 2 0 


cos 0 Q 
-TTT * 


/tan | 


sin^0 s 


/ tan — 

S = -K 2 csc 2 0 - K 2 2£Lfs + log 2 

sin 2 0 s e W a n"£ 


itan^ 
\ 2 , 


sin z 6 sin*0 


dS = 2K2ml- ^A + log 0 fcl 


esc 0 d0 


S dR = SR - \ R dS 


csc‘ 2 0 + 


cos 0 g 
sin 2 0 c 


! + log. 


0\~l 2 l 0 


/ Q\ 

cos 0 cosjls / tan 2 

sin 2 0 sin 2 0 S 6 Wan £§ 


2K esc 0 d0 
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APPENDIX - Continued 


Substituting equations (All) and (A12) into equation (A10) gives the following equa- 
tion for the first integral: 



\ 2K esc 9 dB 

J e w 


_(l - V 2 )sin 6 


2K esc 6 d 6 


Pw e w 


Substituting equation (A13) into equation (A9) gives 
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APPENDIX — Concluded 


Solving the second integral and rearranging gives 


'D,a 


yM 0 


P t,W 7fC 

(l - V 2 )sin 0 

jPoo Al 

V 

L 

_ 


■ V w 2 )sin 6> w 2K y 

— fit lc ®el 


/tan @-\ 

_ ! + fbJ 

i . 1 

\tan jy 

AL J 


Recognizing that 


and 


"w _ _w_ 
A L W oo 


n£ v(l - V 2 ) 7 ' 1 
Aj_, sin 6 


C_ - 2 \J-L 

c D,a o ' — 

rwj 


w 


Poo 

v_ 


r 


(l . v 2 ) 1 " 1 - JL(l . vj )*- 1 


2Ky v(l - V 2 ) 7 " 1 
y - 1 sin 8 


' tan ^ 


log, 


■el. 0 W , 

tan — 


1 + — 

w 


Since 


then 


C p - 


2 

V 

r M oo 2 

_Poo 


iw(l - v 2 ) y ’ 1 - 1 


W ^ l , w v \jl - v / t ^ 

D >a - C » - C P>» - (y - 1)M J P. sin ? 

2 > 


4K p t,w v(l - V 2 ) 7 " 1 


tan- 


(A14) 
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Pressure coefficient , 



Ray angle , 6 , deg 


Figure 2 - Pressure coefficient in conical flow field. IVL = 2; y = 1.4. 
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Figure 6.- Variation of additive-drag coefficient and mass-flow ratio with ray angle. 9 S = 30°; y = 1.4. 
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